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ABSTRACT 

We calculate the expected angular power spectrum of the temperature fluctuations 
in the microwave background radiation (MBR) generated in the quasi-steady state cos¬ 
mology (QSSC). The paper begins with a brief description of how the background is 
produced and thermalized in the QSSC. We then discuss within the framework of a 
simple model, the likely sources of fluctuations in the background due to astrophysical 
and cosmological causes. Power spectrum peaks at Z « 6 — 10, 180 — 220 and 600 — 900 
are shown to be related in this cosmology respectively to curvature effects at the last 
minimum of the scale factor, clusters and groups of galaxies. The effect of clusters is 
shown to be related to their distribution in space as indicated by a toy model of struc¬ 
ture formation in the QSSC. We derive and parameterize the angular power spectrum 
using six parameters related to the sources of temperature fluctuations at three charac¬ 
teristic scales. We are able to obtain a satisfactory fit to the observational band power 
estimates of MBR temperature fluctuation spectrum. Moreover, the values of ‘best fit’ 
parameters are consistent with the range of expected values. 

Subject headings: microwave background inhomogeneities, cosmology 


1. INTRODUCTION 

The quasi-steady state cosmology (QSSC) has been proposed by Hoyle, Burbidge and Narlikar 
(1993, 1994a, b , 1995) as an alternative to the standard hot big bang model. This cosmology does 
away with the initial singularity, and does not have any cosmic epochs when the universe was very 
hot. The synthesis of light nuclei and the origin of the microwave background radiation (MBR) 
have therefore been explained by different physical processes than those invoked in the hot big 
bang. (cf. Hoyle, Burbidge and Narlikar, 2000; Burbidge and Hoyle, 1998). 
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We concentrate here on the MBR and our purpose is to explain its observed fluctuations in 
terms of the large scale features of the QSSC. First we begin with a brief description of how it is 
generated in this cosmology. A typical QSSC scale factor for a flat Robertson-Walker universe is 
given by 


S(t) = e^ p [l + r] cos(27t t/Q)], (1) 

where the time scales P ~ 10 3 Gyr 3> Q ~ 40 — 50 Gyr are considerably greater than the 
Hubble time scale of 10-15 Gyr of standard cosmology. The function r(t ) is very nearly like t, 
with significantly different behaviour for short duration near the minima of the function S(t). 
The parameter r/ has modulus less than unity thus preventing the scale factor from reaching zero. 
Typically, r/ ~ 0.8 — 0.9. Hence there is no spacetime singularity, nor a violation of the law of 
conservation of matter and energy, as happens at the big bang epoch in the standard model. This 
is because, matter in the universe is created through minibangs, through explosive processes in the 
nuclei of existing galaxies that produce matter and an equal quantity of a negative energy scalar 
field C. Such processes take place whenever the energy of the (7-field quantum rises above the 
threshold energy mpc 2 , the restmass energy of the Planck particle which is typically created. For 
details of the basic physics see Hoyle et al (1995) and for models of the kind (1) see Sachs et al 
(1997). 

While the overall energy level of the (7-field is below this threshold, it can rise in the neigh¬ 
bourhood of highly collapsed massive objects (objects close to becoming black holes), and these 
then become the sites of creation. Also since the overall energy density of the C— field goes as S' -4 , 
the creation is facilitated near the epochs of smallest S. It is close to these relatively denser epochs 
that the creation activity is at its peak, switching itself off as the universe continues to expand in 
the cycle. Hence from one cycle to another the density of matter would have fallen by the factor 
exp (-3Q/P), but for the compensating creation of matter at the beginning of a new cycle at the 
minimum of S. 

Thus in the simplest model, although the universe is in a long term steady state, it oscillates 
over shorter time scales, with each cycle physically the same as the preceding one. Each cycle has 
the new matter taking part in the process of formation of stars and galaxies which evolve exactly 
as in the previous cycle. What happens to the radiation of stars in this process? Assuming that 
the radiation density falls like S~ 4 in each cycle, it would be depleted by the factor exp(—4 Q/P) 
from one oscillatory minimum to the next. It is this gap that is made up by the starlight generated 
during the cycle, since the universe is in a steady state from one cycle to the next. Given P and 
Q, and the starlight energy generated per cycle, we can estimate the radiation background at any 
stage of a cycle. We have shown elsewhere (Hoyle, Burbidge and Narlikar 1994a, 2000) that it is 
this starlight thermalized by dust grains which is responsible for the microwave background (MBR) 
in the QSSC. Further, the estimate of starlight generated during the typical cycle allows one to 
show that after thermalization the temperature at the present epoch should be very close to 2.7 K 
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(Burbidge and Hoyle 1998; Hoyle et al 2000). 

In the following sections we first review the spectrum and homogeneity of the MBR in the 
QSSC, and then describe the likely causes of inhonrogeneity of the MBR. We finally compare the 
results with recent observations in which the fluctuations of the MBR at different angular scales 
have been measured. 


2. THE SPECTRUM AND TEMPERATURE OF THE MBR 

In this model it is hydrogen burning in stars in galaxies that generates the energy of the MBR, 
and the energy distribution of stellar radiation over different wavelengths from a typical galaxy 
varies through the cycle in the following way. At the minimum S epoch, new gaseous material is 
acquired, and new stars are formed, with the energy coming out mostly in the blue to ultraviolet, 
from stars of masses greater than M & . However, as the cycle proceeds towards the maximum S 
phase after the elapse of times of the order ~ 25 Gyr, the massive stars burn out, leaving only the 
low mass stars still shining. The typical stars thus left are the dwarfs of types K to M , with the 
consequence that the radiation is then mainly in the red and the infrared. However, as the cycle 
proceeds towards the next minimum, these wavelengths are shortened since the universe contracts 
by a considerable factor. 

The total energy density of starlight generated in a given cycle can be estimated from the 
present observations of the starlight background and by extrapolating the star burning activity over 
the full cycle. Suppose that this energy density is e. If the background radiation energy density u 
at the start of a cycle (when the scale factor was at its minimum value S m ) were u m , then in the 
absence of any new addition to it, following the u oc S~ 4 law, this would drop to it m exp(— A.Q/P) 
at the end of the cycle. To maintain a steady state from one cycle to next therefore the shortfall 
in u m is made up by e. Hence 


e = u m {l - e 4Q/P }. (2) 

This fixes the value of u m in terms of starlight energy input. If we know the maximum redshift 
z m in the present cycle, we find that the present energy density of the MBR is 

u 0 = u m (1 + z m ) . (3) 

The details of e, z m , P and Q were discussed in the papers cited above (see, for example, Hoyle 
et al 2000). A present day MBR temperature of ~ 2.7 K arises from such considerations. However, 
to talk of a ‘temperature’ one must first show that the MBR has reached thermal equilibrium. We 
next describe briefly what must occur if this result is to be achieved. 
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To fix ideas, we will consider the typical oscillatory cycle from one minimum of scale factor to 
next. Let 


x = S(t)/S(t min ) (4) 

where, in going from one minimum to the next we have ignored the exponential term. It is clear 
from this time dependence, that in going from the maximum to the minimum, the scale factor 
drops by a factor (1 + r])/( 1 — rf) ~ 12 for 77 ~ 0.85. We shall assume as an additional parameter 
that the present epoch is characterized by a; = 6 . 

Thus, infrared wavelengths of 1000 nm at x = 6 are shortened to ultraviolet wavelengths of 
only ~ 160 nm at the x = 1 (minimum S ) epoch. However, by this stage, the intergalactic density 
is increased by a factor 6 3 as the minimum is approached and the issue of absorption becomes 
important. For it is through frequent absorptions and re-emissions that the radiation background 
gets thermalized. 

We have shown that metallic particles in the form of graphite whiskers as well as iron whiskers 
play a crucial role in the absorption process. As Narlikar et al (1997) have discussed, a plausible 
case based on laboratory experiments and astrophysical evidence, can be made for the condensation 
of metallic whiskers from the hot ejecta of supernovae, which are blown out into intergalactic space 
by the shock waves generated at the time of explosion. The extinction properties of such whiskers, 
typically of length ~ 1 mm and diameter 0 . 01 /im differ considerably from those of normal spherical 
dust. In particular the iron whiskers at cryogenic temperatures are a dominant source of opacity 
in the microwave region, while carbon whiskers are more effective at the shorter UV wavelengths. 

Typically a galaxy belongs to a cluster and we expect these whisker grains to fill up the 
intergalactic space within the cluster. The starlight from a galaxy in the cluster will therefore 
pass through such a dust as it emerges out of the cluster. We expect therefore absorption and 
re-radiation of starlight outside galaxies and also, in a larger angular scale, outside clusters. As 
we shall see later, the production of microwaves in this fashion will go on in each cycle, and the 
process of frequent absorption and re-radiation by whiskers will eventually generate a uniform 
background, except for the contribution from the latest generation of clusters. These will stand out 
as inhonrogeneities on the overall uniform background. 

The value of absorption coefficient Q a bs for graphite whiskers is essentially constant for all 
wavelengths longer than ~ 1 /irn, extending even to the long radio wavelengths, being equivalent to 
an absorption coefficient of 10 5 cm 2 g 1 . However, it is three times this value for the UV radiation. 
Now an intergalactic density of ~ 1CU 34 g cm ' 3 at x = 6 would rise to ~ 2 x 10~ 32 g cm -3 , at the 
minimum (x = 1) epoch. Over a cosmological distance of 10 27 cm at this stage, the optical depth 
would be ~ 6 for UV, and ~ 1 for wavelengths longer than 1 p, m. 

It is significant that this type of dust plays a role in the observed redshift apparent magnitude 
relation of Type la supernovae. As shown by Banerjee et al (2000), an intergalactic dust density of 
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the above order gives a very good fit to the observations (Perlmutter et al 1999, Riess et al 1998). 
Since the results are very sensitively dependent on dust density, it cannot be fortuitous that the 
density for ‘best fit’ in the supernova case is of the right order to explain thermalization of the 
MBR. Further, recent studies of the high redshift supernova SN 1997// (z = 1.7) show the model 
to be consistent with observations within the error - budget currently applicable (Vishwakarma 
2002, Narlikar et al 2002). 

The great bulk of the optical radiation that becomes subject to thermalization in the contract¬ 
ing phase of the oscillation will have traveled a distance of the order of 10 Gpc, or even more in the 
case of microwaves. The radiation incident on a carbon whisker will have been in transit since the 
maximum S epoch of the previous cycle, and it includes all the microwave radiation existing before 
the present cycle, as well as the starlight generated by the galaxies in the current cycle. The result 
is that all this radiation is well thermalized and uniform in energy density. However, the carbon 
whiskers themselves will be lumpily distributed, on the scale of clusters of galaxies. So, as the 
minimum of S is approached at the end of the last cycle, the conversion of starlight to microwaves 
will have been lumpy. As the starlight is progressively absorbed, the typical grain temperature T g 
first rises above the MBR temperature Tmbr before dropping back to it as the grain radiates. 

The effect of grains being lumpily distributed on the MBR is to cause the slight rise followed by 
the fall back to the MBR temperature to be lumpy too. The radiation background itself, however, 
does not have this lumpiness: because the total assembly of grains has a negligible heat content 
and in thermal equilibrium, each grain emits as much heat as it receives. Since the emissivity of 
the particles does not depend on wavelength, each emits a Planckian spectrum corresponding to 
its temperature T g . Eventually when all particles come down to the temperature TmbRi however, 
further absorption and reemission bring about a strict Planckian distribution at this temperature. 
Although as was suggested by Narlikar, Wickramasinghe and Edmunds (1975), the thermalization 
can be achieved by the carbon whiskers, the presence of a small quantity of iron whiskers helps the 
process further. The absorption and re-radiation by the whiskers at the oscillatory minimum, thus 
generate a mixing of radiation from distances as far as ~ 10 29 cm, at the relatively low intergalactic 
particle densities prevalent there, so in effect permitting radiation to travel freely. 


3. THE ORIGIN OF MBR INHOMOGENEITIES 

This picture suggests that the overall background will be very smooth. However, there will be 
some tiny fluctuations of intensity stamped on it arising from certain intrinsic inhomogeneities of 
the process as well as the cosmological model. We consider them in decreasing order of angular 
scale. 

For this purpose we propose a simple model in which the sources of inhomogeneities can be 
traced to the epoch of last minimum of oscillation (x = 1 in our present case). This is the maximum 
density epoch when the extinction by dust will be most effective. The characteristic geometrical 
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size of the model is determined by the spacetime curvature R and will be of the order R -1 / 2 . Any 
large scale inhomogeneity of the MBR would arise over this characteristic cosmological size. This 
size features in standard cosmology also and is referred to as the ‘horizon size’ or the ‘Hubble radius’ 
(see Weinberg 1972). For smaller sizes we need to take note of inhomogeneity on the scale of rich 
clusters of galaxies, since they represent local concentrations of starlight not yet fully thermalized 
and redistributed. On still smaller scales we expect to see inhomogeneities on the scale of groups 
of galaxies or even individual galaxies. (There should likewise be some inhonrogeneity on the scale 
of superclusters, but it is expected to produce a weaker effect and we will ignore it in the present 
model.) 

It is clear from the above that the issue of inhomogeneities of the MBR in the QSSC is linked 
with the way large scale structures are formed in this cosmology. A structure formation scenario 
in the QSSC has not been developed to the level of sophistication that it has acquired in standard 
cosmology. A toy model by Nayeri et al (1999) has shown, however, that the structure formation 
process as in the QSSC is very different from that in the standard cosmology, being essentially 
driven by the minicreation events with gravitation playing a secondary role. Thus existence of 
galaxies and clusters at high redshifts (z ~ 5) is taken for granted in the QSSC picture of structure 
formation. The toy model of Nayeri et al shows how clustering at various scales develops from an 
initial random distribution, through creation events and expansion. In a ‘steady state’ situation 
the model has to demonstrate how the physical conditions in one QSSC cycle are repeated in the 
next. We will elaborate on this idea in the following section. 

To give quantitative estimates of the above effects, we will choose the same cosmological pa¬ 
rameters that we have used to explain the redshift magnitude relation based on Type la supernovae 
(see Banerjee, et al, 2000). They are : 

Ho = 65 km s -1 Mpc -1 , z max = 5, Ao = — 0.358. (5) 

The parameters P and Q do not enter explicitly in the calculation but we may typically take 
them to be 1000 Gyr and 50 Gyr, respectively. 

The largest angular scale will arise from the size of spacetime curvature. At the present epoch 
this will be of the order of ( c/H 0 )~ 2 . This curvature corresponds to a linear size c/H 0 . As we go 
towards the minimum scale factor epoch, this scale decreases. The angle subtended by the scale at 
us from the minimum scale epoch can be worked out. For the parameters P ~ 1000 Gyr, Q ~ 50 
Gyr, z max = 5 — 6, we get an angular size of ~ 10°. We therefore expect this scale to show up in 
the regime of relatively large angular scale fluctuations of the intensity of MBR. The first discovery 
of inhonrogeneity of MBR by COBE (Smoot et al 1992) was of this order. 

To assess the magnitude of fluctuations, we reproduce the arguments of Hoyle, at al (1994a). 
Consider the whisker grain distribution at the epoch of last minimum of S(t). The whiskers may 
be distributed non-uniformly to begin with, with some regions having larger whisker densities than 
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others. Until the deviations from uniformity of radiation background become too small, they are 
able to push the grains down the temperature gradients to restore uniformity. Consider a fluctuation 
by a small factor y in the average energy density ~ 5 x 1CU 10 erg cm -3 of the microwaves at this 
stage. A fluctuation of this order is able to move grains of density p g and velocity V provided 
1/2 p g V 2 = 5 x 10 -1 V With p g ~ 2 x 10 -33 g cm -3 and V ~ 0.1c (an adequate speed to fill a 
region of size ~ 10 26 cm in the available cosmological times of the order 10 17 s), we get y ~ 2 x 10~ 5 , 
corresponding to AT/T ~ 5 x 10~ 6 . This is of right order as found by COBE. Thus any temperature 
fluctuations above this value would create temperature gradients of such magnitude as to move the 
grains around from high to low density regions and to thereby create a uniformity. Only small 
enough fluctuations would therefore remain intact over regions of size rsj 10 26 cm. Regions of this 
order are characteristic of the curvature size referred to above. 

The strongest fluctuations will, however, come from rich clusters of galaxies lying in the present 
cycle. These denote the late additions to the background already in existence from starlight of 
previous cycles. Note that the light from stars belonging to galaxies in clusters of earlier generations 
will already have been fully thermalized and such clusters will not stand up as fluctuations. It is 
only the newer clusters at the last oscillatory minimum with inadequate time to have been merged 
in the rest of the thermal background, that would produce fluctuations. A typical cluster produces 
the extra starlight, which on subsequent thermalization produces additional temperature AT over 
the average background temperature T. Such a fluctuation will be confined to the neighbourhood 
of the cluster. We show next how to estimate it. 

Imagine a 10 Mpc size region, typically containing ~ 10 4 galaxies, each emitting starlight at 
the rate ~ 10 44 erg s~ 4 (corresponding to an average galaxy with absolute magnitude —21.4). The 
flux of radiation in the form of degraded starlight across the surface of this region (assumed to be 
a sphere of radius 5 Mpc), will be ~ 3.35 x 10 -4 erg s _1 cm -2 . Now, the present energy density of 
MBR ps 4.2 x 10 13 erg cm 3 , implies that by the (1 + z) 4 rule, it was 6 4 times this value at the 
epoch of the last minimum of S. Using Stefan’s law the flux across a sphere of radius 5 Mpc for 
this radiation will be c/4 times the energy density, i.e., about 4.08 erg s _1 cm" 2 . Comparing the 
excess flux due to the galaxies in this region computed above with this value we find that there is 
an excess flux equal to 8.21 x 1CU 5 of the average background. Equating this to 4A T/T, we find 
that the temperature fluctuation AT = 56pK. 

Because of the inhomogeneity of the universe on the scale of clusters, the calculation cannot be 
made more precise, and at this stage one should only look at its order of magnitude. The fact that 
this calculation yields a temperature fluctuation of the order reported by the various observations 
is, from the point of view of the QSSC, very encouraging. 

Clusters of galaxies are normally associated with a temperature decrement because of the 
Sunyaev - Zel’dovich effect. Indeed observations show such effects for clusters of redshifts < 1. By 
contrast, the effect described above is predominant for clusters of large redshifts ~ 5, corresponding 
to the epoch of the last minimum of the scale factor. Thus we expect the thermalization effect 



to dominate over the S-Z effect for such clusters, whereas for later epochs the former effect gets 
reduced since dust extinction is much less. A more detailed study of cluster formation and evolution 
in the QSSC is required to decide at what stage clusters begin to acquire hot gas and when the gas 
temperature rises to levels when the S-Z effect becomes important. After such a study is carried 
out it will be possible to estimate the average S-Z decrement of temperature in any given direction 
by integrating over redshifts going up to the above stage. 

Over and above the clusters we expect smaller scale fluctuations on the scales of individual 
galaxies. The metallic whiskers are generated and ejected from inside a typical galaxy. They have 
escape - velocities that take them well beyond the galaxy into its immediate external environment. 
Thus we expect that the extra starlight generated by the galaxy will provide a slight excess of AT 
in this region after thermalization. The typical length scale of a small group of galaxies would be 
1 - 2 Mpc. 


4. RELATIONSHIP TO STRUCTURE FORMATION 

From the above discussion it is clear that according to the QSSC the fluctuations of the MBR 
arose at the relatively recent epochs (associated with the last minimum of expansion of the scale 
factor) and they would be related to the large scale structure in the universe. We will therefore 
outline briefly the current ideas on this topic in the framework of the QSSC. We will consider in 
particular the formation and distribution of clusters of galaxies, since they will turn out to have 
the most significant effect on the smoothness of MBR. For details we refer the reader to Banerjee 
and Narlikar (1997) and Nayeri, et al (1999). 


4.1. Gravitational Stability 

To begin with, it is necessary to contrast the structure formation process in the QSSC with that 
in the standard cosmology. In standard cosmology, structure formation begins in the form of small 
scale fluctuations in the spacetime metric as well as matter contents, fluctuations which are believed 
to be of quantum origin. These fluctuations then evolve as they grow under the effect of gravitation, 
with the inflationary phase playing an essential role by modifying the spectrum of inhomogeneities 
to the scale-invariant form. The inhomogeneities continue to grow with gravitational clustering 
playing a major role. Since radiation and matter (at least the baryonic part) were strongly coupled 
in the early epochs, the growth of fluctuations in matter affects radiation also. This interaction 
continues till the surface of last scattering and the inhomogeneities imprinted on the radiation 
background at that stage get imprinted on the MBR and can be observed today. 

It is essential to appreciate that the above scenario does not occur in the QSSC whose past 
history is quite different from that of the standard cosmology. In the QSSC the quantum gravity 
dominated epoch as well as inflation did not occur; nor was there a surface of last scattering. Thus 
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to understand the presence of fluctuations in the MBR one needs to look at an entirely different 
scenario in which creation of matter at periodic intervals of a finitely oscillating universe having a 
long-term de Sitter type expansion, plays a key role. 

That the gravitational force does not play the primary role in structure formation in this 
scenario is seen from the work of Banerjee and Narlikar (1997), wherein these authors looked at the 
evolution of small departures from the spacetime metric of the QSSC as well as its matter density 
and flow vector. Specifying the unperturbed metric by the flat Robertson-Walker line element, 

ds 2 = dt 2 - 5 2 (t)[(dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 ] (6) 

with S(t) given by (1), the perturbed metric is written as 

9n u = T hfn,'), = ho^, <700 = 1 T /loo- (7) 

Here rj is the Minkowski metric (p, v = 1,2, 3; x M spacelike coordinates and x° = t) and | |<C 

1. Likewise the density is perturbed from po to po + pi, \ pi/p |<C 1 and the flow vector from 
Uq = (1,0, 0,0) to (uq + u\) with | u\ |<C 1. 

The perturbed set of field equations therefore describe the full gravitational effect on the 
perturbations. Banerjee and Narlikar (op.cit.) found that in a typical oscillation,these perturbations 
grow to a limited extent before subsiding. These authors therefore concluded that gravitational 
effects cannot play a major role in forming large scale structure in the QSSC. 

In fact the dynamics of the QSSC requires the existence of an additional force (besides grav¬ 
itation) which is derived from a scalar field C (Hoyle et al, 1995, Sachs et al, 1996). This scalar 
field in turn is related to the creation of matter. A look at the field equations of the QSSC will help 
understand why this field acts in such a way that gravitation does not play a key role in structure 
formation. These equations are: 


Rik 2 9 ik^d T A(7 ik ki 


M (c) 

Tik + 


M _ (C) 

where T is the usual matter tensor while T A is given by 


( 8 ) 


(C) 

T ik = -f 


CiC k 


\gikC l C e 


( 9 ) 


Both k = 87 rG/c 4 and / are positive. Thus the energy tensor of the C-field has a negative coefficient 
(—/) which makes its effect on spacetime structure repulsive, i.e., acting in direction opposite to 
gravity. Note that this tensor acts even when no creation is going on. This is why, during a typical 
oscillation the normal gravitational processes of clustering are impeded. 
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A parallel may be cited from standard cosmology with a positive cosmological constant (A > 0). 
In those models which have the A-term sufficiently strong even in the early stages, gravitational 
clustering and formation of structures cannot take place. For example, with the presently favoured 
value of A, clustering ceased at z ~ 1. 

However, the C-field intervenes in a different way through its property that it induces creation 
of matter. As discussed next, matter is created near collapsed massive objects (i.e., objects close to 
becoming black holes) and the created matter is ejected outwards because of the negative stresses 
generated by the C-field. It is this creation phenomenon that holds the key to structure formation 
in the QSSC. 


4.2. Matter Creation : A Toy Model 

The field equations of the C—field show that at the point of creation of a particle of momentum 
Pi, the gradient of the C—field satisfies the condition C,; = pi, (i = 0,1, 2,3). Thus, if the rest mass 
of the particle created is mo, then (with the speed of light c = 1) we have at creation C;C* = m q. 

The use of this relation in the Schwarzschild solution near a massive object of mass M shows 
that at a coordinate distance r the magnitude of the C-field energy is raised to 


C,C n 


o 

m 

Y 2GM ’ 
c 2 r 


( 10 ) 


where m 2 is the magnitude of the energy at large distances from the massive object. Now it may 
happen in general that m < mo, so that no creation takes place in general. However, near a highly 
collapsed massive object for r close enough to the Schwarzschild radius of the object, the creation 
condition is satisfied and creation can take place (Hoyle and Narlikar 1966). 

The process normally begins by the creation of the C-field along with matter in the neighbour¬ 
hood of a compact massive object. The former, being propagated by the wave equation, tends to 
travel outwards with the speed of light, leaving the created mass behind. However, as the created 
mass grows, its gravitational redshift begins to assert itself, and it traps the C-field in the vicinity 
of the object. As the strength of the C-field grows, its repulsive effect begins to manifest itself, 
thus making the object less and less bound and unstable. Finally, a stage may come when a part 
of the object is ejected from it with tremendous energy. It is thus possible for a parent compact 
mass to eject a bound unit outwards. This unit may act as a centre of creation in its own right. 

We shall refer to such pockets of creation as minibangs or mini-creation events (MCEs). A 
spherical (Schwarzschild type) compact matter distribution will lead to a spherically symmetric 
explosion whereas an axi-symmetric (Kerr type) distribution would lead to jet-like ejection along 
the symmetric axis. Because of the conservation of angular momentum of a collapsing object, it 
is expected that the latter situation will in general be more likely. Using this picture Nayeri et al 
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(1999) proposed a toy model for structure formation and evolution, which is summarized next. 

A large number of points (N ~ 10 5 — 10 6 ), each one representing a mini-creation event, is 
distributed randomly over a unit cubical area . The average nearest-neighbour distance for such a 
distribution will then be L = N 1//3 . Now suppose that in a typical mini-creation event, each point 
generates another neighbour point at random within a distance, d = xL in 3D . Here, the number 
x is a fraction between 0 and 1 . We shall call x the separation parameter. It denotes an ejected 
piece lying at a distance < d from the parent compact object. 

The sample cube is then homologously stretched by a linear factor 2 1 / 3 to represent expansion 
of space. We now have the same density of points as before, i.e., 2N points over a volume of 2 
units. From this enlarged cube remove the periphery so as to retain only the inner unit cube.This 
process thus brings us back to the original state but with a different distribution of an average N 
points over a unit cube. This process is repeated n times . Here the number of iterations, n, plays 
the role of “time” as in the standard models of structure formation. The number distribution of 
points evolves as the ‘creation + expansion process’ generates new points near the existing ones. 

Not surprisingly, soon after, i.e., after n = 3 — 4 iterations of the above procedure, clusters and 
voids begin to emerge in the sample volume and create a Persian Carpet type of patterns. As the 
experiment is repeated, voids grow in size while clusters become denser. Fig. 1 illustrates the point 
distribution within a typical thin slice of thickness 0.001 of the cube and shows that expansion 
coupled with creation of matter is a natural means of generating voids and clusters. 

To bring the toy model closer to the reality of the QSSC, Nayeri, et al proceeded as follows. 
Since the creation activity is expected to be confined largely to a narrow era around a typical 
oscillatory minimum, when the C-field is at its strongest, by considering the number density of 
collapsed massive objects at one oscillatory minimum of QSSC to be /, the number density at 
the next oscillatory minimum would fall to / exp (—3 Q/P), that is, if no new massive objects were 
added. Thus to restore a steady state from one cycle to the next, within each unit volume 


uf=[l- exp (—3 Q/P)\f ~ (3 Q/P)f, (11) 

masses must be created anew. In other words, a fraction u of the total number of massive objects 
must replicate themselves in the above fashion. 

Notice that, unlike the old steady state theory which had new matter appearing continuously, 
we have here discrete creation, confined to epochs of minimum of scale factor. The ‘steady-state’ 
is maintained from one cycle to next. Which is why the above fractional addition oj is required at 
the beginning of each cycle. 

Therefore, instead of creating a new neighbour point around each and every one of the original 
set of N points, one does so only around u)N of these points chosen randomly, where the fraction 
u is as defined in equation (11). Likewise, the sample volume is homologously expanded by the 
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N= 1,000,000 x = 0.8 n= 10 Az = 0.001 



X Direction 


Fig. 1.— A simulation of large-scale structure in the QSSC, based on the creation of one generation 
of clusters in the vicinity of those of the previous generation, keeping aligned ejection from one 
generation to the next. Filaments and voids develop on going from one generation to next. 
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factor exp (3 Q/P) only instead of by factor 2 . 


5. The Two Point Correlation Functions 

Although visual inspection of Figures like 1 suggests that the toy model is proceeding along the 
right lines, a quantitative measure of the cluster-void distribution helps in comparing simulations 
with reality. The dimensionless autocorrelation function 

£( r ) =< [ p { r )- < p >][p(ri + r)- < p >] > / < p > 2 , (12) 

where < p > is the average density in the volume, is one convenient measure of such irregularities 
in the space distribution. Typically, different classes of objects cluster at different characteristic 
lengths. To fix ideas Nayeri et al (1999) looked at distribution of clusters of galaxies. Observation- 
ally, it is believed that the two point correlation function for cluster distribution obeys the following 
scaling law: 


&c(r) = (^) > ( 13 ) 

with 7 ~ 1.8 and ro = 25/i _1 Mpc, where the Hubble constant at the present epoch is taken to 
be 100/j kms _ 1 Mpc _1 . In order to quantify the issues of formation of structures in this scenario 
Nayeri, et al (op.cit) took the following measures. 

It is known that instead of having a uniform distribution of matter on large scales, the observed 
universe has structures of typical sizes of a few tens of megaparsecs . These “structures” are regions 
of density considerably higher than the background density, with the maximum density contrast 
5 = ( 5p(r)/ < p >) going from order unity (in the case of clusters) to a few thousand (in the case 
of the galaxies). 

Any process which generates structures must be able to produce to the zeroth order, entities 
whose density contrast is of such magnitude and with the property that on larger and larger distance 
scales, the density contrast became less significant. This ensures that on a large enough scale the 
universe is homogeneous. 

Given this prescription for generating structures without gravitational dynamics, Nayeri, et al 
first ensured that the visual impression created by the cluster simulations did imply that as the 
number of iterations were increased the number of high density regions also increased. In the initial 
(totally random) configuration one expects to find regions of high density arising only because of 
the Poisson noise. In the later “epochs” after a few iterations, however, one would expect to find, 
as in a clustering scenario, that the variation of the one point distribution function for density 
(p/ < p >) with < p > the average density in the volume, showed a steady and significant increase 
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lter=7 

lter=5 



-2 -1.5 -1 -0.5 0 

log(r) 


Fig. 2.— A two-point correlation function for simulations of the kind shown in Fig. 1 has a power 
law-type distribution with index approaching —1.8 after a few iterations of the process of creation 
The evolution of the index towards the ideal —1.8 slope line (shown by the continuous line) is show] 
as the process goes on for up to ten generations. 
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in the the number of high and intermediate density regions. This expectation was borne out and it 
was also observed that the value of maximum density also increased as a function of the number of 
iterations, which in this experiment corresponds to “time”. The density field was generated on a 
grid placed into the simulation volume using the algorithm of cloud in cell. The simulations showed 
the growth of structures through rise in the density maximum as a function of number of iterations. 

A quantitative measure that was computed from these data sets was the two-point correlation 
function. Fig. 2 shows the two point correlation function for the case of the QSSC based model. 
As “time” goes on, the slope of the correlation function gets closer and closer to —1.8. From the 
value of the X-axis intercept of the two point correlation function Nayeri, et al could get a rough 
estimate of the size of the structures in units of the size of our simulation box. From their results 
they estimated that the size of the structures formed is approximately (3 = 0.15 — 0.3 times the 
boxsize. If one sets these values equal to the observationally accepted value of ?’o, one can get a 
better physical sense of the results. If we set, /3 = 0.3, say, and ?’o = 25 h~ l Mpc, then the linear 
size of the simulation box would be ~ 84/r -1 Mpc. Typical cluster masses and sizes resulted from 
such a numerical exercise. 

This may be seen as an attempt to relate the toy model to a realistic cosmological scenario. Of 
course, as the above exercise shows, the results can be scaled up/down by rescaling the simulation 
parameters and thus are independent of the ‘absolute’ size of the box. A more detailed dynamical 
theory of the creation process will help relate the absolute size of clustering to the theoretical 
parameters. 

We now examine the typical angular scales of these three types of inhomogeneities and how 
they are reflected in the angular power spectrum of the MBR. So far as clusters are concerned, we 
will fold in the above evolving two-point correlation function into our calculation. 


6. THE ANGULAR POWER SPECTRUM OF INHOMOGENEITIES 

Typically we will look at the angle subtended by a linear scale a of an object located at redshift 
2 i in the quasi-steady state model whose scale factor is defined in equation (1). To compute the 
above angle, we define the dimensionless parameters by the following formulae: 


Ao 

fI c o 

K 0 


8nGpo 

~3Hf 


density parameter, 


ml 

8nGp c o 

ml 


cosmological constant parameter, 


creation density parameter, 


H?,S{ 


curvature parameter. 


(14) 
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Here A is the cosmological constant, which is negative in the QSSC. In the model considered by 
Banerjee et al (2000) as well as here, Kq = 0. The angle subtended at the observer by the above 
scale at redshift z (see Banerjee and Narlikar, 1999) is then given by 


a = 


H 0 a( 1 + z) 


rl+z 


dy 


n -1 


' i \/Ao + y 3 + H c0 y 4 


(15) 


We will take z = z max = 5. (The maximum redshift used here is indicative only. The results do 
not change much even if it were increased to 6-10.) The calculation, using the parameters of (5), 
then gives 


7ra 

a =-. 

3600 

where a is measured in units of megaparsecs. Unless explicitly stated otherwise, angles will be 
referred to in radians. 

This then is the characteristic angle subtended by an object of linear size a Mpc, viewed when close 
to the last minimum. Assuming that this is the size of a rich cluster of galaxies, we expect that 
the characteristic angular size of the MBR inhomogeneity it generates will be of this order. As we 
will shortly show, this does translate to a peak in the angular power spectrum at a harmonic of 
order l ~ 1/a. So we arrive at the following approximate relationship between cluster size and the 
harmonic it generates: 



l 


3600 

ira 


(17) 


If we set a Mpc ~ 2x radius of curvature at the last minimum, we expect a large angle peak at 
l ~ 6 — 10. The smaller angle peaks will occur at much larger values of l. Thus a peak in the power 
spectrum of MBR anisotropy at harmonic l ~ 200 as observed by the Boomerang and the Maxima 
groups (de Bernardis et al 2000, Hanany et al, 2000) is generated in this theory by rich clusters of 
diameter ~ 5.5 Mpc, located at cosmological distances. Smaller peaks are expected on the scale of 
groups of galaxies, at the higher values of l « 500 — 1000, and on the scale of superclusters at the 
lower values Z ~ 10 — 20. These latter peaks will be smaller than the cluster peak at l ~ 200 and 
thus they will be much harder to detect. 

We next show that these expectations are borne out by detailed calculations of this effect on 
the angular power spectrum of MBR. We model the above effect along the lines of Hajian and 
Souradeep ( 2002 ) as follows. Imagine a typical type of inhomogeneity as a set of small disc-shaped 
spots, randomly distributed on a unit sphere. The spots may be either ‘top hat’ type or ‘Gaussian’ 
type. In the former case they have sharp boundaries whereas in the latter case they taper outwards. 
We will assume the former for clusters, and the latter for the curvature effect and also for galaxies 
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or groups of galaxies. This is because the clusters will tend to have rather sharp boundaries whereas 
in the other two cases such sharp limits do not exist. 

Let us begin with the simple case of a random distribution of infinitesimally small spots. With 
each direction hi for a typical spot centre on the sky chosen randomly, we have from N such spots 
of roughly identical temperature fluctuation A T c 

N 

A T{fi) = ^2 ® n i3(h - hi) A T c . (18) 

Here h is a general direction on the unit sphere. If the locations of the spots in the sky are 
uncorrelated, we have 


(fl ni Ojuj) — $ij- ( 19 ) 

Writing AT as a function on the sphere, we expand it in a series of Yi m (8, (f>): 


AT(h) = 22 a im Y im(n). 

lm 


( 20 ) 


The expansion of 5{h 


hi) in terms of the spherical harmonic functions allows us to obtain 


N 

a lm. = y ] ■ ( 21 ) 

2=1 

The corresponding angular power spectrum is given by 


_ 1(1 + 1 ) 

= 2t r 


<21 + 1 ) a lm) 

v J m 


1(1 + 1) N 

2ir 47r 


(A T c ) 2 . 


( 22 ) 


In the last step we use eq. (19) that encodes our simplistic assumption of uncorrelated spots. We 
get, for large l , Ci oc l 2 . This result is valid for angular scales much larger than the angular 
scale of the typical spot. The l 2 rise is a direct consequence of assuming that the location of the 
spots are uncorrelated. More detailed modeling where the correlation between the spots is included 
would serve to soften the l 2 rise. It may be an important issue when comparing our results with 
observations of MBR inhomogeneity on large angular scales. 

However, we need to model spots of finite sizes ranging from ~ 10° to a few arc minutes, 
corresponding to the three typical inhomogeneities discussed above. For spots of finite size, assumed 
to be circularly symmetric with a given profile centered on the typical spot centre, we have a 
modification of (18) in the form (where the delta function is replaced by a function /): 
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N 

AT(n) = ^o„AT c /(n • h/). (23) 

8=1 

Using the circular symmetry of the spots, we expand 

fifi -hi) = ^ 21 ^ fi Pi(h • hi), (24) 

i 

in terms of its Legendre transform, f). It is then possible to express 

c i = > fi Pint ■ hj). (25) 

i j 

Assuming the location of the N spots to be uncorrelated and randomly distributed (see eq. (19)) 
this gives 


W 2 2 

1 2vr 4 tt v cJ Jl 


(26) 


For a normalized Gaussian spot profile with typical (l-<7 width) angular size a , the function 
fl is given by 


fi(a) = exp (-l 2 a 2 / 2). (27) 

The angular power spectrum Ci oc l 2 e~ l ~ a2 has a single peak at l p « a -1 . 

For a top hat profile, wherein the spot has uniform temperature fluctuation across a finite disc 
of angular size a, with a sharp drop-off to the background temperature outside, the function fi is 
given by 


fi{a) 


1 

ww 


cos a Pi (cos a) — Pi -1 (cos a) 
(2 sin 2 a/2) 


(28) 


The Ci in this case has a series of peaks in the multipole space. The first peak in C; for top-hat 
spots occur at l p ~ it/ 2a and has the largest amplitude. The successive peaks occur roughly at 
integer multiples with diminishing peak amplitude. 


We therefore consider a composite model in which the inhomogeneities of the MBR arise from 
a superposition of random spots of three characteristic sizes corresponding to the three effects 
discussed above. From the order of magnitude estimates made earlier, we expect three peaks in the 
angular power spectrum at their respective l - ranges. To check the accuracy of this expectation, 
we consider a 6 - parameter modeling of the angular power spectrum 
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Ci 


Ai l(l + l)e- l 2 a '+A 2T ^ [ 

+A 3 1(1 + l)e- /2 “i 


cos a 2 Pi(cos a 2 ) — P;_i(cos a 2 ) 2 
2 sin 2 0 : 2/2 


(29) 


with the parameters {A\, A 2 , A 3 , aq, a 2 , 03 } determined by obtaining the best fit with the present 
observations of the angular power spectrum of MBR inhomogeneities. The parameters A \, A 2 and 
A 3 depend on the number density as well as the typical temperature fluctuation of each kind of 
spot. However, given Ai, A 2 and A 3 , it is possible to compute the rms temperature fluctuation 
contributed by each component to the overall MBR temperature fluctuations. 

If the model is successful, the values of a\,a 2 ,a 3 (or, the corresponding multipole value l p at 
which the C; from each component peaks,) and the rms temperature fluctuation for components 
would turn out to be within the range of values expected on the basis of the physics of the processes. 

The analysis till now has assumed that the distribution of hot spots is uncorrelated. As 
described in the previous section, toy model of structure formation within QSSC does recover the 
observed power law spatial correlation in the distribution of clusters. We extend the computation 
of CMB anisotropy spectrum to the case when the hot spots are correlated. We do this only for 
the hot spots that are linked to the clusters. Then eq. (19) is revised to 


(«ni«n j) = W(0) (30) 

where w{9) is angular correlation of clusters on angular scale, 6 = cos _1 (nj.nj). The form of Ci in 
eq. (26) is revised to 

A N 

Cl = -^l(l + l)/ f V (31) 

where w(9) = ^(2Z + l)/(47r)'U;P;(nj.n ? -). The angular correlation w{9) can be related to the 3-D 
spatial correlation function £(r) using the well known Limber equation ( see Peebles, 1980). In 
particular, for a power law correlation function £(r) = (r/r 0 ) -7 , it can be shown that w(0 ) oc 9 l ~^ 
for 9 <C 1 limit and 1 < 7 < 6. In this case, U[ oc for 7 < 3 (Peebles and Hauser, 1974). 
The correlation of the hot spots due to clusters are expected to be bounded within 7 = 3 , the 
limit of uncorrelated spots and 7 ~ 2, the correlation function that the clusters distribution evolves 
to within the QSSC model. Consequently, we also consider an extended model with an extra 
parameter 7, corresponding to an angular power spectrum 


Ci 


Hi 1(1 + l)e- /2 “? + A 2 ^ 

+A 3 l{l + l)e- /2 “i 


cos o^PHcos a 2 ) — Pi- i(cos a 2 ) 2 
2 sin 2 07/2 


( 32 ) 


where all the other parameters are as described in eq. (29). 
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Fig. 3.- The band power spectrum of the MBR temperature inhonrogeneity, measured (> 2a 
detections only) by the different experiments. The compilation of Podariu et al. 2001 has been 
used. 
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7 . COMPARISON WITH OBSERVATIONS 

We now compare our ‘three component model’ with a published data compilation available 
at the time of writing. Fig. 3 describes the error-bars of the observed data which indicates the 
extent of scatter at present (Podariu et al. 2001). The data have been binned into sixteen bins in 
multipole space with a mean AX) = \/Ci over the bin and a 1 — a error-bar (Podariu et al. 2001, 
table 2). Two peaks are easily visible at l ~ 10 and 200 while a third peak at the much higher value 
of ~ 600 is probably present. The l - values can be related to physical dimensions of the sources of 
inhomogeneities by the formula (IT). For the actual fitting, we have used the above binned data. 

We first discuss the case of uncorrelated spots (7 = 3). Fig. 4 shows the ‘best-fitting’ angular 
power spectrum curve AT; obtained for QSSC by using a six parameter model for the characteristic 
angular sizes, a* of the spots and the corresponding amplitudes, A % . For a given type of spot, the 
angular size can be readily related to the multipole l p where the contribution to C; peaks. The 
amplitudes are used to compute the rms temperature fluctuation, AT rms , contributed by each of 
the components. For the ‘best fit’ model we find l p \ = 27, l p 2 = 182, l p 3 = 502, AX rms i = 20 /jA, 
AT rmS 2 = 31/cA and AT rmS 3 = 16/cA. For the QSSC model we have used six parameters to obtain 
the above fit with a minimum y 2 of 16.8. We divide the minimum y 2 value by 16 — 6 = 10 degrees 
of freedom to arrive at the value 1.68 for the reduced y 2 for our fit. This is useful for a comparison 
with any other model. 

It is encouraging to note that the parameters for the best fit model are broadly consistent with 
our expectations. The value of l p 3 = 500 is somewhat below the expected range of 600 - 900. This 
is most likely because the largest central multipole value of the binned data set is l ~ 600. As MBR 
anisotropy on large multipole bands are determined and included in the binned data, we expect the 
best fit values of l p 3 (and possibly, AT rmS 3 ) to go up. We also note the significance of the top-hat 
approximation to spot profile in clusters in our fit to the data. If we replace the top-hat profile of 
cluster anisotropy by a Gaussian then we can find a best fit model with y 2 ~ 23 with parameters 
l p i = 25, l P 2 = 201, l P 3 = 900, AT rms i = 19/cA, AT rms 2 = 30/rAT and AT rmS 3 = 20/cA. Thus the 
reduced y 2 ~ 2.3. The similarity in the amplitude and sizes clearly demonstrates the robustness 
of our model to finer details of the spot profile. The reality may lie between the two limits, the 
Gaussian profile and the top-hat one. 

We compared with the same binned data, the anisotropy spectrum prediction of a grid of 
open-CDM and A-CDM models within the standard big bang cosmology. We varied the matter 
density, flo = 0.1 to 1 in steps of 0.1; the baryon density, Q^h 2 from 0.005 to 0.03 in steps of 0.004 
where h is the Hubble constant in units of 100 km s _1 Mpc -1 ; and the age of the model, to from 
10 Gyr to 20 Gyr in steps of 2 Gyr (Sugiyama, 1995). For each value of flo we considered an open 
model and one where a compensating was added to get a flat model. For the same binned data 
set, we found the minimum value of y 2 = 23 for the flat universe model with Ho = 0-7 (Ha = 0.3), 
1 \h 2 = 0.021, to = 14 Gyr (h = 0.52). Recent MBR inhomogeneity data, when combined with 
high redshift supernovae data that support non-zero and constraints on Hq from large scale 
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Fig. 4.— Best-fitting angular power spectrum curves in the QSSC for 7 = 3 (continuous curve) 
and 7 = 2 (dotted curve) are plotted and are compared with the one for the favoured big bang 
model with 17 a = 0.7 (dashed curve). The fitted data set is the one described in Figure 3 and has 
been averaged into 16 bins in multipole space (Podariu et al. 2001). 
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structures observations, favour flat cosmological models with higher values of I2 a- In our analysis, 
these models have higher, but comparable, y 2 values, eg., for VL^h 2 = 0.021 and to = 14 Gyr, the 
models with 12 a = 0.6 (h = 0.62) and 12 a = 0.7 (h = 0.67) have y 2 = 24 and 25, respectively. For 
a comparison, the corresponding curve for the favoured big bang model (12a = 0.7) is also shown 
in Fig. 4. 

Next we incorporate the correlation in the spots arising from clusters corresponding to spatial 
correlation in the distribution of clusters. In the absence of a detailed modeling of CMB anisotropy 
arising from clusters, we admit the entire range 3 > 7 > 1.8. For 7 ~ 2, we find the best fit values 
l p i = 18, l P 2 = 230, l P 3 = 977, AT rms i = 13/r/i, AT rmS 2 = 32/iK and AT rmS 3 = 23/j.K. For these 
best fit parameters, y 2 = 39. As seen in Fig. 4, most of the discrepancy arises at the observed peak 
around l ~ 200, particularly from the 13th point whose contribution to y 2 alone is ~ 37 percent. 
This point also does not fit the standard big bang models and seems to be a general outlier. With 
7 = 2 the second peak in our model is broad and cannot fit the high amplitude of the peak observed 
in the binned compilation. At present, this may not be a significant problem for our model even 
with 7 = 2 . As 7 is increased to 3, the fit improves rapidly. We are not in a position to fix the 
value of 7 that would be appropriate for the clusters at z ~ 5 within QSSC. On the observational 
front, the amplitude at the peak needs to be well established by an experiment such as MAP that 
covers both the low l and high l parts with good Z-space resolution. The plot of the individual 
band power estimates in Fig. 3, shows considerable dispersion in the central value of band power 
estimates and some of discrepancies are attributed to calibration uncertainty. 

In Fig. 5, we compare the band power values \ZC~i, averaged within the 16 multipole bins 
for the best fit QSSC and the favoured big bang (12 a = 0.7) models, with the binned values of 
the observed MBR inhomogeneity data. Within QSSC, the unclustered hot-spots allow for a very 
good fit to the CMB data up to l 600. The peak around l ~ 200 in the present data is best fit 
by uncorrelated spots and fit around the peak is adversely affected as the correlation between the 
spots is increased. We also find that the details of spot profile does not affect the fit significantly. 


8. CONCLUDING REMARKS 

In the QSSC the MBR arises from the matter in galaxies and this is quite different from its 
origin in the hot big bang cosmology. For example, the major peak (at l ~ 200) in the MBR in the 
framework of the QSSC is explained in terms of rich clusters of galaxies, whereas in the hot big 
bang cosmology it is the ‘Doppler’ peak associated with acoustic oscillations of the photon -baryon 
fluid. The QSSC interpretation links the inhomogeneities of the radiation field to those of the 
matter field, visible matter in the form of clusters of galaxies. Indeed, one may argue that in the 
QSSC interpretation, the observations of MBR-inhomogeneities provide us with a direct diagnostic 
of the structural hierarchy in the universe. For example, the pattern of patchiness around l ~ 200 
depends on the cluster - cluster correlation function. However, we must use the clusters at large 
redshifts (~ 5) for accurately estimating the magnitude of the effect. In this paper we account for 
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effect of cluster-cluster correlation function based on earlier published results obtained using a ‘toy 
model’ structure formation model within QSSC. We hope to carry out such an analysis when more 
complete data on the cluster distributions is available going back to epochs of high redshifts and 
the theory of structure formation and evolution in the QSSC makes further progress. 

Additionally, from Fig. 4 we see that peaks occur at different locations for l > 500, for the most 
favoured big bang model, and the best-fitting QSSC model. Future observations at high values of 
l from MAP, PLANCK and other experiments may help distinguish between the two cosmologies. 

In the end perhaps it is best to stress the attitudinal difference between the two cosmologies. 
In the big bang cosmology, the inferences are related to the postulated initial conditions prevailing 
well beyond the range of direct observations (at redshifts > 1000). In the QSSC the attempt is to 
relate patchiness of structure, (at redshifts « 5) which may be observable one day, to the patchiness 
of MBR. These latter studies admittedly do not give predictions as sharp as those given by the 
former, but they may perhaps claim to be less speculative. 
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